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ABSTRACT 

Recent studies suggest that only three of the twelve brightest satellites of the Milky 
Way (MW) inhabit dark matter hales with maximum circular velocity, Knax, exceeding 
~ 30 km/s. This is in apparent contradiction with the ACDM simulations of the 
Aquarius Project, which suggest that MW-sized halos should have at least 8 subhalos 
with yniax > 30 km/s. The absence of luminous satellites in such massive subhalos is 
thus puzzling and may present a challenge to the ACDM paradigm. We note, however, 
that the number of massive subhalos depends sensitively on the (poorly-known) virial 
mass of the Milky Way, and that their scarcity makes estimates of their abundance 
from a small simulation set like Aquarius uncertain. We use the Millennium Simulation 
series and the invariance of the scaled subhalo velocity function (i.e., the number of 
subhalos as a function of v, the ratio of subhalo V^iax to host halo virial velocity, 
V200) to secure improved estimates of the abundance of rare massive subsystems. In 
the range 0.1 < ly < 0.5, A'sub(> ^) is approximately Poisson-distributed about an 
average given by (Nsuh) ~ 10.2 (j//0.15)~^'^^. This is slightly lower than in Aquarius 
halos, but consistent with recent results from the Phoenix Project. The probability 
that a ACDM halo has 3 or fewer subhalos with T^nax above some threshold value, 
Vth, is then straightforward to compute. It decreases steeply both with decreasing Vth 
and with increasing halo mass. For Vth = 30 km/s, ~ 40% of Mhaio = 10^^ Mq halos 
pass the test; fewer than '^ 5% do so for Afhaio^2 x 10^^ Mq; and the probability 
effectively vanishes for Mhaio^3 x 10^"^ Mq. Rather than a failure of ACDM, the 
absence of massive subhalos might simply indicate that the Milky Way is less massive 
than is commonly thought. 
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1 INTRODUCTION 

The striking difference between the relatively flat faint-end 
slope of the galaxy stellar mass function and the much 
steeper cold dark matter halo mass function is usually rec- 
onciled by assuming that the efficiency of galaxy forma- 
tion drops sharply w ith decreasing halo mass (see, e.g., 
IWhite fc Frenklll99lh . Semi-analytic models of galaxy for- 
mation have used this result to explain the relatively small 
number of luminous satellites in the Milky Way (MW) halo, 
where ACDM simulations predict the existence of thou- 
sands of subhalos massive enough, in principle, to host dwarf 
galaxies. In these models, the small number of MW satel- 
lites reflects the relatively small number of subhalos massive 
enou g h to host luminous ga l axies (see, e.g. , 
1 19931 : iBullock et ai] l200(3l : iBenson et aU 



20021: ICooper et al.ll2010l:lLi et al.ll2010l : iMaccio et allboiol : 
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Quo et al.ll2011l : lFont et al.ll201ll ). 

This is a model prediction that can be readily tested 
observationally, given the availability of radial velocity mea- 
surements for hundreds of stars in the dwarf spheroidal satel- 
lites of the Milky Way. Combined with photometric data, 
radial velocities tightly constrain the total mass enclosed 
within the luminous ra dius of these satellites (jWalker et al.l 
I2OO9I : IWolf et aLlbOld ). The latter correlates strongly with 
the total dark mass of the dwarf, which is usually expressed 
in terms of its maximum circular velocity, Vmax, a quantity 
less a ffected than mass by tidal stripping (|Penarrubia et al.l 
120081 ). 

Kinematical analyses of the Milky Way dwarf 
spheroidals have been attempted by several authors in re- 
cent years, with broad consensus on the results, at least 
for the best-studied nine brightest dwarf spheroidal MW 
companions: Draco, Ursa Minor, Fornax, Sculptor, Ca- 
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rina , Leo I, Leo II, Canis Venatici I, and Sextans (see , 
e.g..|Pena rmbia ct al."2008i : IStrigari et al.ll2008l: lLokas"2009l: 
IWalker et al., i2009i: iWolf et all i20ld : IStrigari et al.i i201of). 
These studies suggest that some of these galajcies may in- 
habit halos with V^ax as low as 12 km/s, and agree that 
allj appear to inhabit halos with values of Vmax below a low 
threshold, Kh ~ 30 km/s. Only three dwarf irregular satel- 
lites - the Magellanic Clouds and the Sagittarius dwarf - 
may, in principle, inhabit halos exceeding this threshold. 

The most straightforward interpretation of this result 
is that massive subhalos in t he Milky Way are rare . How- 
ever, as argued recently by iBovlan-Kolchin et al.l (|201ll . 
120121 '). this is at odds with the results of the Aquarius 
Project, a series of N-body simulations of six different ha- 
los of viriaEI mass in t he ran ge 0.8 < Ivhoo/IQ^'^ Mq < 1.8. 
iBovlan-Kolchin et al.l (|201ll ) noted that the largest subha- 
los in these simulations are significantly denser than inferred 
for the halos that host the brightest dwarf spheroidals in the 
Milky Way. 

As discussed by iParrv et al.l ((201^) and 
iBovlan-Kolchin et al.l (|2012l l. the discrepancy can be traced 
to the fact that the largest Aquarius subhalos are signifi- 
cantly more massive or, equivalently, have too large a value 
of Vinax to be compatible with the measured kinematics of 
the brightest dwarf spheroidals. Specifically, the Aquarius 
halos have, on average, ~ 8 subhalos with Vinax > 30 km/s 
within the virial radius, larger than the Knax of the bright- 
est dwarf spheroidals, prompting questions about why these 
massive subhalos fail to host luminous satellites in the Milky 
Way. If this result holds, it may point to a failure of our ba- 
sic understanding of how galaxies populate low mass halos 
or, more worryingly, of the ACDM paradigm itself. 

Two issues may affect these conclusions. One is that 
the Aquarius Project simulation set contains only 6 halos 
and, therefore, estimates of the abundance of rare massive 
subhalos are subject to substantial uncertainty. The second 
point is that the number of massive subhalos is expected to 
depend sensitively on the virial mass of the host halo, which 
is only known to within a factor of 2-3 for the Milky Way. 

We address these issues here by using large numbers 
of well-res olved halos identified in the Millennium Simula - 
tion series (jSpringel et al.l2005l : lBovlan-Kolchin et al.ll2009r ). 
This is possible because, in agreement with earlier work, we 
find that the abundance of subhalos, when s caled appro- 
priately, is independent of halo mass (see^ e.g.,lMoore et al.l 



10^1 



ll999l : lKravtsov et al.ll2004l : ISpringel et allboOSl ). We use this 
to derive improved estimates of the average number of mas- 
sive subhalos, as well as its statistical distribution. The prob- 
ability that a halo has as few massive subhalos as the Milky 
Way can then be evaluated, both as a function of host halo 
mass and/or subhalo mass threshold. 

This paper is organized as follows. Sec. [2] describes 
briefly the simulations we used in our analysis. We present 



^ One possible exeeption is Draeo, where the data might allow a 
more massive halo. 

^ Unless otherwise noted, we define virial quantities as those cor- 
responding to spheres that enclose a mean overdensity A = 200 
times the critical density for closure. M20O1 for example, corre- 
sponds to the mass within the virial radius, T2Q0- When other 
values of A are assumed the subscript is adjusted accordingly. 
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Figure 1. The number of subhaloes with Vmax ^ Vth as a func- 
tion of the virial mass of their host haloes, M200, in the Millen- 
nium Simulations (MS and MS-II), as well as in the levcl-2 runs 
of the Aquarius and Phoenix Projects. Subhalos are identified 
within the virial radius, r200i of their host systems. Difi'erent sym- 
bols correspond to each simulation, as labelled, and are coloured 
according to the value of the threshold, Vt\-i- Error bars denote 
the rms plus Poisson error in each mass bin. Note the nearly lin- 
ear dependence of the number of subhalos with halo mass. Due 
to numerical resolution, few subhalos with velocities less than 
~ 100 km/s are found in the MS simulation, so the Vth = 30 
and 60 km/s MS curves in this case are omitted for clarity. For 
massive, well-resolved halos the results are much less affected by 
numerical limitations and there is good agreement between MS 
and MS-II. Subhalo abundance is insensitive to small variations 
in the cosmological parameters. Triangles connected by a dotted 
line show results corresp onding to a run that adopted the lat- 
est WMAP7 parameters I IKomatsu et al.ll201ir) ; in contrast, the 
Millennium Simulations adopted parameters consistent with the 
Ist-year analysis of WMAP data. 



our main results in Sec. [3l and end with a brief summary in 
Sec.H 



2 SIMULATIONS 



The two Mi Uennium simulations (MS: [Springel et allbOOSl 



(,MSj 
l2009l ) 



and MS-II; iBovlan-Kolchin et al.l 120091 ) provide the main 
datasets used in this study. Both are simulations of a flat 
WMAP-1 ACDM cosmogony with the following parameters: 
^M = 0.25, Q.b = 0.045, h = 0.73, n, = 1 and ag = 0.9. 

The MS run evolved a box 500 Mpc//i on a side, with 
2160^ particles of mass m^ = 8.6 x IQ^Mq/H. MS-II evolved 
the same total number of particles in a box 1/125 the volume 
of MS and had, therefore, 125 better mass resolution {rrip = 
6.885 X IQ^Mq/K). The nominal spatial resolution is given 
by the Plummer-equivalent gravitational softening, which is 
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Figure 2. The number of subhalos with maxiinuirL circular ve- 
locity exceeding a certain fraction, u, of their host halo virial 
velocity, as a function of N200 > the number of particles within the 
virial radius of the host. Curves for three values oi u = 0.1, 0.15, 
and 0.2 are shown. Error bars for MS and MS-II indicate the 10 
and 90 percentile in each bin, and are omitted when the bin con- 
tains a single halo. The excellent agreement between results for 
MS and MS-II at given Af200 reflects the halo mass invariance of 
the A'aub(> 1^) function; each particle is 125 X more massive in 
MS than in MS-II. Results converge for well-resolved halos (i.e., 
those with large N20o)- As expected, the smaller u the larger the 
minimum number of particles, Afjgg", needed to obtain converged 
results. 



Figure 3. The scaled subhalo velocity function, i.e., the number 
of subhalos with maximum circular velocity exceeding a certain 
fraction, u = Vniax/V200i of the host halo virial velocity. Dotted 
and dot-dashed curves show averages for the six Aquarius halos 
and nine Phoenix halos, respectively. Dashed and solid curves 
correspond to MS and MS-II. Four curves are shown for each, 
corresponding to averages over all halos in mass bins of width 
0.1 dex centred at log^g M200/MQ = 11.5, 12.5, 13.5, and 14.5. 
Error bars (shown only for the lowest and highest mass bins) 
indicate the rms scatter in each bin. Only halos satisfying the 
constraint N200 > Af^g"(v) are used. All simulations arc in good 
agreement when well-resolved halos are considered. The scaled 
subhalo velocity function is thus nearly invariant with mass. See 
text for further discussion. 



ep = 5 kpc/h, and 1 kpc//i for the MS and MS-II runs, 
respectively. 

We also use halos from the Aquariu s Project 
llSpringel et al.ll200a ) and the Phoenix Project (|Gao et al.l 
12012! ) (level-2 resolution). These are ultra high- resolution 
simulations of six MW-sized halos (M200 ~ 10^^ Mq) and 
nine cluster-sized halos (M200 ~ W^^Mq), each resolved 
with a few hundred million particles within the virial ra- 
dius. 

The normalization of the power spectrum adopted in 
these simulations is slightly high er than favoured by th e lat- 
est WMAP dataset (WMAP7: [Xomatsu etlj] I2OIII ). but 
this is expected to affect the abundance of halos of given 
virial mass rather than the mass function of subhalos, which 
is the main focus of our study. We have verified this explic- 
itly by analyzing a 1620^-particle simulation of a 70.4 Mpc/h 
box that adopts the WMAP7 cosmological parameters (see 
Fig.[T}. The particle mass in this run is 6.20 x 10^ Mq//i and 
gravitational interactions were softened with ep = 1 kpc/h. 

Halos and subhalos are id entified in all simulations by 
SUBFIND (|Springel et al.|[200ll ). a recursive algorithm that 
identifies self-bound structures and substructures in N-body 
simulations. 



3 RESULTS 

We first investigate the scale invariance and other statistical 
properties of the distribution of subhalo Vmax and then apply 
our results to subhalos in the Milky Way. 



3.1 Subhalo Knax distribution 

Fig. [1] shows, as a function of host halo virial mass, the to- 
tal numbeilfl of subhalos with maximum circular velocity, 
Vmax, exceeding a specified velocity threshold, Vth- Results 
are shown for three different values of Kh- The average num- 
ber of subhalos in each halo mass bin is shown by symbols 
connected by solid (MS-II) or dashed (MS) lines. Individual 
level-2 Phoenix and Aquarius halos are shown by crosses 
and open squares, respectively. WMAP7 results are shown 
by open triangles connected by a dotted line. 

Fig. [1] illustrates that: (i) the number of subha- 
los depends roughly linearly on halo mass and increases 
strongly with decreasing velocity threshold, and that (ii) 



Unless otherwise noted, we identify subhalos within the virial 
radius, r^oo, of the host halo. 
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the slight change in cosmological parameters from WMAPl 
to WMAP7 has a negligible effect on subhalo abundance. 

Fig. [T] also shows that numerical resolution limits the 
halo mass and velocity threshold for which convergence in 
subhalo abundances is achieved. Indeed, there are fewer sub- 
halos in MS, the simulation with poorest mass resolution; so 
few with velocities less than ~ 100 km/s that the Vth = 30 
and 60 km/s MS curves have been omitted for clarity. When 
halos and subhalos are resolved with enough particles, how- 
ever, the results converge well. For Kh = 120 km/s, MS, 
MS-II, and WMAP7 halos yield similar numbers of subhalos 
over the whole halo mass range considered, despite the fact 
that, at given M200, MS-II and WMAP7 halos have ~ 125 x 
more particles than their MS counterparts. Furthermore, the 
results for Phoenix and Aquarius are in good agreement with 
MS-II, even though halos in MS-II have 700 x fewer parti- 
cles than Aquarius and 2x fewer particles than Phoenix, 
respectively. 

We explore the requirements for numerical convergence 
in more detail in Fig. [2l where we plot, as a function of the 
total number of particles within the virial radius, A^200, the 
average number of subhalos with Vmax exceeding a certain 
fraction of the host halo virial velocity: A'aub(> v), for p = 
Knax/V2oo = 0.10, 0.15, and 0.20. Results are shown for MS 
and MS-II halos with dashed and solid curves, respectively. 

This figure highlights two important points. One is that 
at given p there is good agreement between all simulations 
provided that halos are resolved with enough particles. The 
second is that, when the first condition is met, A'sub(> p) is 
independent of halo mass. (Recall that, at fixed A^200, MS 
halos are 125 x more massive than their MS-II counterparts.) 
This agreement, together with the fact that the A^sub(> p) 
curves plateau at large values of N200, iniply that the scaled 
subhalo velocity function (i.e., the number of subhalos as a 
function of !/ = Knax/V2oo) is invariant over many decades 
in halo mass. Fig. [5] also makes clear that numerical conver- 
gence requires that a halo be resolved with a total number of 
particles above some ((/-dependent) minimum number, N^q 
(listed in Table 1). The converged values agree well, within 
the statistical uncertainty, with the results for Phoenix and 
Aquarius halos. 

We show the invariance of subhalo abundance explic- 
itly in Fig. [21 where Afsub(> i') is plotted for MS and 
MS-II halos grouped in 4 bins of different halo mass. 
Only halos satisfying the A''2oo > N^^^ constraint are 
used here. T hese res ults confirm earlier work (see , e.g. , 
Moore et al][l9 99: K ravtsov et al.ll2004l: [Zheng et al.ll2005l : 
Weinberg et al.ll2008l : ISpringel et al.ll2008l '). and imply that 
we can combine all well-resolved halos into one large sample 
to derive robust estimates of the statistical distribution of 

A'"aub(> U). 

This is shown in Fig. [4] for p = 0.1, 0.15, and 0.2, com- 
puted using aU MS and MS-II halos with iV2oo > A''20o"('^). 
as given in Table [T] In the top panel, which corresponds 
to subhalos identified within r2oo, the histograms show the 
A'aub(> y) distributions for the 614, 3070, and 6867 ha- 
los that satisfy, respectively, the minimum particle num- 
ber constraint. The bottom panel shows subhalo numbers 
identified within a slightly larger radius, rioo, which is on 
average ~ 30% larger than r2oo- (The mean and rms dis- 
persion of the distributions of A^sub(> p) are listed in Ta- 
ble 1.) Note that the results obtained for MS and MS-II 
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Figure 4. The distribution of N^^^,{> p) for p = 0.1, 0.15, and 
0.2, computed for well-resolved MS and MS-II halos; i.e., those 
with particle numbers exceeding N^^{p) (as given in Table [l]|. 
The top panel refers to all subhalos within the virial radius, r20o; 
the bottom panel to subhalos within a radius roughly 30% larger, 
rioo- The average and rms for Aquarius and Phoenix halos are 
shown at the top of the plot. Note that subhalos in Aquarius seem 
slightly overabundant relative to either Phoenix or the Millen- 
nium Simulations, but still well within the statistical uncertainty. 
The Afsub(> p) distribution is well approximated by a Poisson 
distribution: the solid curves show Poisson distributions with the 
same averages as each histogram. 



are in excellent agreement with the results obtained for 
Phoenix halos. Subhalos in Aquarius are slightly overabun- 
dant, perhaps because of small biases in their assembly his- 
tories (JBovlan-Kolchin et al. 201(i), but still consistent with 
the MS and MS-II results given the large variance {a% ) 
expected from the sample of only 6 Aquarius halos. Note as 
well that, as expected, the larger volume encompassed by 
rioo yields larger subhalo numbers than found when identi- 
fying subhalos only within r2oo. For r < r^oo, the average 
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Figure 5. Probability that a halo contains 3 or fewer subhalos 
with Vmax > 30 km/s, as a function of halo mass (top tickmarks) 
or virial velocity (bottom tickmarks). The solid black curve cor- 
responds to assuming Poisson statistics and that A^sub(> '^) = 
10.2(i//0.15)~^'^^ , the average number of subhalos within the 
virial radius of well-resolved MS and MS-II halos with particle 
numbers exceeding A'^^q" (see Table [ij . The sensitivity of the re- 
sult to the assumed minimum number of particles is shown by the 
red and blue curves, which correspond to increasing the values of 
A^^q"(i/) by factors of 4 or 8, respectively. Results using only the 
nine Phoenix or six Aquarius halos are shown in cyan and green, 
respectively. Note that because subhalos are slightly overabun- 
dant in Aquarius (see Fig. |4]l the probabilities arc systematically 
lower than when considering either Phoenix halos or the Mil- 
lennium simulations. The same is true if subhalos are identified 
within a radius larger than the virial radius. The dashed curve 
shows probabilities when the search radius around each halo is 
increased by roughly 30% to rioo- 



-'Vsub(> 1^) is a steep function of v, well approximated, in 
range 0.1 < i^ < 0.5, by 



the 



(iVsub)(>i') = 10.2(i//0.15)' 



-3.11 



(1) 



Fig. |4] also shows that the distribution of A'aub(> v) fol- 
lows Poisson statistics closely; the solid curves are actually 
not fits, but just Poisson distributions with the same av- 
erage as each of the histograms. Clearly, these provide a 
good description of the distribution of A*'sub(> v) at fixed 
u. This conclusion i s supporte d by earlier work (see, e.g., 
iKravtsov et al.ll2004l : iBovlan-K olchin et al.lboiOl '). as weU as 
by the data listed in Table [l] the average number of sub- 
halos is roughly similar to the variance, as expected from a 
Poisson process. 



3.2 Massive satellites in the Milky Way 

We can use these results to address the Milky Way missing 
massive satellites problem highlighted in Sec. [T] In particu- 



lar, it is straightforward to compute the probability that a 
halo has X or fewer subhalos with Knax > 30 km/s within 
its virial radius, once a virial mass (or, equivalently, a virial 
velocity, V200) has been assumed for the Milky Way. This is 
given by. 



X 



fc=0 



(2) 



where A^ = (A'sub)(> v) is given by Eqn.[T] 

The solid black curve in Fig. [S] shows /(^ 3) as a func- 
tion of virial mass (upper tickmarks on the abscissa) or virial 
velocity (lower tickmarks). The probability is a steep func- 
tion of the assumed halo mass: more that 40% of W^^Mq ha- 
los pass this test, but only ~ 5% of 2xlO'^^M0 systems do so. 
The probability becomes negligible for M2oo^3 x 10^^ Mq. 
This suggests that the scarcity of massive subhalos is best 
thought of as placing a strong upper limit on the virial mass 
of the Milky Way, rather than as a failure of the ACDM 
scenario. 

It is important to assess the sensitivity of this conclu- 
sion to the parameters assumed in this study. For example, 
should the velocity threshold be placed at 25 km/s , rathe r 
than at 30 km/s, as argued bv lBovlan-Kolchin et al.l l|2012h . 
the upper limit on the mass of the Milky Way would become 
even more restrictive. The results, however, could still be 
read from Fig. [5] after shifting the tickmarks by 30/25 = 1.2 
in the velocity axis or by 1.2^ = 1.73 in the mass axis. 
Thus, for Kh = 25 km/s, a probability of more than 5% 
requires a halo mass M200 < 1 x 10^^ Mq, rather than the 
M200 < 2 X 10^^ M0 appropriate to Kh = 30 km/s. 

We have also examined the dependence of our results on 
-^200" (^) J the assumed minimum number of particles needed 
for convergence (listed in Table [iJ . This is shown by the red 
and blue curves in Fig. \5i which correspond to increasing 






by a factor of 4 and 8, respectively, before deriving 



(A'^Bub)(> 1^)- Fig. [5] makes clear that our results are quite 
insensitive to such changes in A^loff . 

Since Phoenix halos have subhalo abundances in good 
agreement with those in Eqn. [T] our results would not 
change had we chosen the nine Phoenix halos to compute 
{Nsuh){> y) (see cyan curve in Fig. [5]). On the other hand, 
had we chosen to derive (A'sub)(> v) solely from the six 
Aquarius halos, the slight overabundance of subhalos in 
these systems would lead to stricter upper limits on the 
Milky Way halo mass, as shown by the green curve in Fig. [S] 
This result, together with the fact that the average Aquar- 
ius halo mass (1.42 x 10^^ Af©) is uncomfortably close to 
the upper limit discussed above is apparently the reason 
why iBovlan-Kolchin et al.l l|2011f ) originally found such a 
strong discrepancy between the Aquarius simulations and 
the Milky Way. 

Finally, we need to consider the dependence of the num- 
ber of subhalos on the maximum radius used to identify sub- 
structure. The results discussed above refer to subhalos iden- 
tified within the virial radius, r2oo, which is ~ 200 kpc for 
a M200 = 10^^ Mq halo. This is smaller than the maximum 
distance commonly adopted to identify dwarf galaxies as 
Milky Way satellites; for example, Leo I is at roughly 250 kpc 
from the centre of the Galaxy. Therefore, it might be argued 
that subhalos should be counted within a larger radius in or- 
der to make a meaningful comparison. As shown in the bot- 
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torn panel of Fig.|4j subhalos are roughly ^ 50% more abun- 
dant within rioo than within 7-200- For a M200 = 10^^ Mq 
halo, rioo ~ 270 kpc, comparable to the Galactocentric dis- 
tance of Leo I. In analogy with Eqn. 1, the average number 
of subhalos, A'sub(> v), within rioo is well approximated, in 
the range 0.1 < i^ < 0.5, by 



(iVsub)(>!^) = 15.03(1^/0.15)" 



(3) 



The dashed line in Fig.[5]shows that the probability of host- 
ing at most 3 massive subhalos drops significantly when the 
rioo radius is used; only about 20% of M200 = lO^'^ Mq halos 
pass the test then. This stricter constraint emphasizes the 
difficulty of resolving the missing massive satellite problem 
if the Milky Way mass significantly exceeds lO'^^ Mq . 



4 SUMMARY 

We have used the Millennium Simulation series, together 
with the ultra-high resolution simulations of small halo sam- 
ples from the Aquarius and Phoenix projects to study the 
abundance of rare, massive subhalos in ACDM halos. As in 
earlier work, we find that the scaled subhalo velocity func- 
tion (i.e., the number of subhalos as a function of the ratio 
between subhalo maximum circular velocity and host halo 
virial velocity, v = Vniax/V2oo) is independent of halo mass. 
This implies that we can obtain robust estimates of the sta- 
tistical distribution of massive subhalos from large samples 
of well-resolved halos selected from the Millennium simula- 
tions. 

Our main result is that, in the range 0.1 < v < 0.5, the 
number of subhalos within the virial radius, r2oo, is Poisson- 
distributed around an average given by Eqn.[T] Compared to 
this average, subhalos in the Aquarius Project are slightly 
overabundant but still consistent given the large variance 
and the small sample of halos included in that simulation 
suite. Subhalos in the cluster-sized Phoenix Project halos 
are in excellent agreement with Eqn. [l] 

We have then used this result to compute the probabil- 
ity that a halo of virial velocity V200 has a certain number 
of massive subhalos with Vmax exceeding a velocity thresh- 
old, Vth- Applied to the Milky Way, where observations sug- 
gest that no more than 3 (or at most 4) subhalos with 
Knax > 30 km/s host luminous satellites, we find that this 
constraint effectively translates into a strong upper limit on 
the Milky Way halo mass. The probability that a halo with 
Af200^ 3x lO^'^ Mq satisfies this constraint within radius r2oo 
is vanishingly small, but it increases rapidly with decreas- 
ing virial mass. Roughly 45% of M200 = 10^^ Mq halos pass 
this test, and ~ 90% of all halos with Af200 ~ 5 x lO" Mq 
are consistent with the data. These fractions are reduced 
to ~ 20% and ~ 70%, respectively, if subhalos are consid- 
ered within a larger search radius, rioo ~ 1.3r2oo (which, 
for halos of mass ~ 10^^ Mq , is close to the Galactocentric 
distance of Leo I, the most distant bright satellite known 
in the Milky Way). In this case, the number of subhalos, 
(-Waub)(> v), within rioo is given by Eqn. [3] and a Milky 
Way halo mass of M200 = 2 x lO^'^ Mq is strongly ruled out 
by the satellite data. 

The "missing massi ve satellites problem" i n the Milky 
Way hal o highlighted bvlBovlan-Kolchin et all (|201ll . 12013 ) 
and by IParrv et al.l (|2012l l may thus be resolved if the 



mass of the Milky Way halo is ~ 10^^ Mq (see also 
IVera-Ciro et al.l |2012| ). This is well within the range of 
halo masses allowed b y the latest estima tes based on either 
the timing argumen t (ILi fc White! 20081 ) or on abundance- 
matching methods (jGuo et al.l I2OICI ). It is in even better 
agreement with the lower virial masses reported by estimates 
based on (i) the radial velocity dispersion of Milky Way 
satell ites and halo stars (jBattaglia et al.l l2005l : I Sales et al.l 
|2007|): fii) the es cape speed in the solar neighbourhood 
( Smith et al.ll2007l ): or (iii) the kinem atics of halo blue hori- 
zontal branch stars (|Xue et al.ll2008l '). Invoking a ~ 10^^ Mq 
mass for the Milky Way is a simpler and more straight- 
forward resolution than several alternatives advanced in 
recent papers, such as co nsidering the baryon adiabatic 
contraction and feedback (|di Cintio et al.l l2011f ) , reducing 
the central density of subhalos through t idal stripping 
JDi Cintio et all |2012| : IVera-Ciro et all 120121 '). or positing 
radical revisions to the nature of dark matter (jLovell et al.l 
l2012l : IVogelsberger et al]|2012l '). 



We conclude that there is no compelling requirement 
to revise the ACDM paradigm based on the abundance of 
massive subhalos in the Milky Way. There are still, however, 
some uncomfortable corollaries to this solution. One is that 
a 10^^ Mq halo has a virial velocity of only ^ 150 km/s, well 
below the rotation speed of the Milky Way disk, usually as- 
sume d to be Vrot = 220 km/s, or even higher (jReid et al.l 
120091 ). This seems at odds with results from some semi- 
analytic models of galaxy formation that attempt simulta- 
neously to match the TuUy-Fisher relation and the galaxy 
stellar mass function: agree ment with observation seems to 
requi re Kot ~ V200 (see, e.g.. lCole et al.ll2000l : ICroton et al.l 
I2006D . 



A further worry is that an A'f2oo = 10^^ Mq halo might 
not be massive enough to host satellites as massive as the 
Magellanic Clouds. Assuming that Knax for the LMC and 
SMC can be identified with the rotation speed of their HI 
disks, or 60 and 50 km /s, respectively (jKim et al.l 119981 : 
IStanimirovic et al.ll2004 ). we find, using the data in TabIe[Tl 
that only ~ 62% of V200 = 150 km/s halos would be ex- 
pected not to host an LMC-like system. The probability of 
hosting two (or more) subhalos more massive than the SMC 
is of order 20%. None of these probabilities seem unlikely 
enough to cause worry. 



Although our results may explain why few massive sub- 
halos might be expected in the Milky Way halo, this expla- 
nation still assigns MW satellites to very low mass halos, i.e., 
those with Knax < 30 km/s. These halos have masses below 
IO^^Mq, the mass scale below which semi-analytic models 
predict that gal axy formation effi ciency should become ex- 
ceedingly small (|Guo et al.ll2010l ). Given the large number 
of low mass halos expected in a ACDM universe, populat- 
ing even a small fraction of Vmax < 30 km/s systems with 
galaxies as bright as Fornax might lead to substantially over- 
pred icting the number o f dwarfs in the local Universe (see, 
e.g., iFerrero et al.ll201ll ). Without a full accounting of how 
dwarf galaxies form in low-mass halos, concerns about the 
viability of ACDM on small scales will be hard to dispel. 
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Table 1. N2qq(i^) is the minimum number of partieles within the virial radius of a halo needed to achieve convcrgcnee in the abundance 
of subhalos. Af^aios is the number of hales that satisfy such condition in the Millennium Simulations. (N^ub) and "^JV^^,, are the average 
number of subhalos exceeding u and its dispersion, respectively. 



V 

(Knax/Vaoo) 


0.1 


0.15 


0.2 


0.25 


0.3 


0.35 


0.4 


0.45 


0.5 


lymin 
A^halos 


1.0 X lO'' 
614 


2.5 X 10^ 
3070 


1.2 X 10^ 
6867 


7.5 X 10" 
12138 


2.5 X 10-* 
38550 


1.8 X lO"* 
54568 


1.0 X 10"* 
90200 


7.5 X 10^ 
113585 


5.0 X 10-' 
151663 


(A^sub) (r < r20o) 

'^N^ub (^ < ''200) 


36.55 
8.92 


10.14 
3.87 


4.20 
2.27 


2.12 
1.54 


1.14 
1.10 


0.71 
0.85 


0.48 
0.68 


0.34 
0.57 


0.25 
0.48 


(A^subXr <noo) 

'^N^uh (^ < ''100) 


51.91 
12.2 


15.22 
5.23 


6.3 
2.99 


3.11 
1.98 


1.73 
1.39 


1.06 
1.08 


0.69 

0.84 


0.46 
0.68 


0.33 
0.58 
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